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ABSTRACT 

The temperature profile of hot gas in galaxies and galaxy clusters is largely determined by 
the depth of the total gravitational potential and thereby by the dark matter (DM) distribution. 
We use high-resolution hydrodynamical simulations of galaxy formation to derive a surprisingly 
simple relation between the gas temperature and DM properties. We show that this relation 
holds not just for galaxy clusters but also for equilibrated and relaxed galaxies at radii beyond 
the central stellar-dominated region of typically a few kpc. It is then clarified how a measurement 
of the temperature and density of the hot gas component can lead to an indirect measurement 
of the DM velocity anisotropy in galaxies. We also study the temperature relation for galaxy 
clusters in the presence of self-regulated, recurrent active galactic nuclei (AGN), and demonstrate 
that this temperature relation even holds outside the inner region of ~ 30 kpc in clusters with 
an active AGN. 

Subject headings: 



1. Introduction 

Dark matter (DM) dominated cosmological 
structures are a direct outcome of hierarchi- 
cal structure formation with a subdominant 
baryon fraction. The baryons can either cool 
fast to form stars or end up as a hot viri 



alized gas (lllees fe Ostrikerl Il977t ISiikl Il977 

" 10051) . 



Birnboim fc Dekell l2003t iKeres et alj|2005h . This 



hot gas is cust omarily obser ved in galaxies and 
galaxy clusters (|Sarazinlll986l ). and can be heated 
through adiabatic compression, shock heating or 
non-gravitational processes and cooled through 
radiative processes. 

The ability of the hot gas to shock heat to the 



virial temperature opens the possibility that one 
can combine the equations governing the dynam- 
ics of the gas and the DM, namely the equation 
of hydrostatic equilibrium and the Jeans equation. 
The simultaneous solution to these equations can 
in principle allow us to determine the DM ve- 
locity dispersion anisotropy, which holds informa- 
tion about the fundamental difference in the way 
DM and baryons equilibrate in cosmological struc- 
tures. In fact, such a method has already been ap- 
plied to galaxy clusters, where numerical simula- 
tions of cluster formation and evolution have been 
used to confirm that the gas equilibrium temper- 
ature is determined thr ough the averaged velocit; 



dispersion of the DM ([Hansen fc Piffarettil 1200 



ltv 
)7; 
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Host et al. l2009h . The resolution of these simu- 
lations, however, did not allow us to probe radii 
smaller than a few hundred kpc. The tempera- 
ture of the cooling gas remains to be determined 
at smaller radii and in the presence of an active 
AGN. Furthermore, it remains unknown if the gas 
temperature in galaxies will obey a similar simple 
relation. 

In this paper we conduct a dedicated compar- 
ison between the gas and DM in a range of sim- 
ulated cosmological objects. The structure of the 
paper is as follows: first we discuss the relationship 
between the gas temperature and the DM veloc- 
ity dispersion. In section [3] we present the results 
of numerical simulations of galaxy formation, and 
confirm that the gas temperature in galaxies ex- 
hibits a bimodal distribution — the hot gas resides 
at the DM temperature and the cold gas cools be- 
low 10 4 K and forms stars. In section@]we present 
the results of numerical simulations of AGN out- 
flows, which demonstrate that even when the gas 
is heated episodically by a self-regulated AGN, it 
still moves towards the DM temperature, which 
it succeeds in reaching already beyond 30 kpc. 
Finally, in section [5] we explain how our findings 
open up for the possibility of measuring the DM 
velocity anisotropy in galaxies, and in section[6]we 
briefly offer our conclusions. 

2. The temperature relation 

It has been known for years that baryons in a 
DM-dominated potential will heat up to a tem- 
perature that is determi ned by the properties 
of the DM mass profile (iRees fe Ostrikerl Il977t 



of t , and the velocity anisotropy (3 — 1 — of/of , 



1981 



Sarazinl 



Cavaliere fc Fusco-Femianol [1978L 
19861: ICavaliere et al.ll2009h R. This is easily seen 
when considering the similarity between the e qua- 
tion of hydrostatic eq uilibrium ( Sarazinll 19861) an d 
the Jeans equation (jBinnev fc Tremaind 1 19871 ) . 
The former depends on the gas temperature, T, 
and through its equation of state on the density, 
/9 gas , while the latter depends on the DM density, 
Pbm, its radial and tangential velocity dispersions, 



1 The corresponding stellar velocity dispersion will also 
be of the same magnitude as the DM velocity disper- 
sion. However, there are no theoretical reason why the 
anisotro pies of the galaxies should be the same as that of 
the DM l lLokas fc Mam on 2003; Woit ak et alll2009f) . 
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where /snip is the averaged mass of the gas parti- 
cles. For a given total mass profile, the hydrostatic 
equilibrium equation, eq. ([1]), has a wide range of 
possible solutions: basically any radial gas tem- 
perature profile is allowed (in principle), since a 
gas density profile can always be constructed to 
fulfill the equation. Furthermore, there is a priori 
no theoretical connection between the gas temper- 
ature and the dark matter dispersions, except that 
the particles are sitting in the same gravitational 
potential. One can, however, use the dark matter 
dispersions to define a "DM temperature" by 



DMft-B _ 1 / 2 



fim p 



(3) 



where a$ and cr^ are the DM tangential disper- 
sion velocities. We emphasize that there is no 
fundamental reason why the temperature of the 
gas, T, should be identical to the DM tempera- 
ture (in praxis the DM only sets an upper limit 
to the gas temperature). We will, however, con- 
sider the ratio between the gas temperature, and 
the temperature derived from the DM dispersions, 
namely 



/■c(r) 



Tk f 



(4) 



and we will discuss a possible universality of this 
ratio. There is, to our knowledge, no simple dy- 
namical argument why this ratio should be close 
to unity, or why it should have a universal radial 
behaviour. For instance, in the adiabatic simula- 



tions presented in lFaltenbacher et al.l (|2007|) it was 
found that k decreases with radius in the inner 
region. We will here address this question using 
galaxy formation including cooling, heating and 
feedback. 

Naturally, gravity does not care about the 
masses of the individual (collisionless) particles, 
which is why their velocity dispersions are equal, 
and not their temperatures. The standard defini- 
tion of the temperature is written in terms of the 
average kinetic energy of particles and describes a 
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system in thermodynamic equilibrium. This ther- 
modynamic limit is not achieved for most particle 
systems with long-range forces such as gravita- 
tional structures, for which this steady state is nor- 
mally not described by simple distribution func- 
tions, e.g., Maxwell-Boltzmann statistics. Thus, 
the DM is not in thermal equilibrium, and techni- 
cally it does not have a "temperature", however, 
the use of this word should not cause any confu- 
sion. 

An important point is that the DM tempera- 
ture, as defined in eq. ([3]), is not only a function 
of the mass profile or potential of the DM, but 
also depends on the dynamical state of the DM 
through its velocity anisotropy, /3(r). The velocity 
anisotropy is much harder to observationally de- 
termine than the mass profile. Hence for a nonva- 
nishing velocity dispersion anisotropy in the DM, 
the DM temperature is not identical to the veloc- 
ity dispersion of the DM along any random direc- 
tion, such as the radial direction, but only to their 
spatial average. 

We will now proceed to analyse the ratio be- 
tween the local gas and DM temperatures defined 
in eq.([3]), and use numerical simulations of galaxy 
clusters or galaxies to test if k is of order unity. 
Consideri ng the results of two different n umeri- 
cal codes (|Kav et al.ll2007t IValdarninill2006l) simu- 



lating both DM and baryonic physics, Host et al 
(|2009l ) demonstrated that the « = 1 within 20% 
for relaxed galaxy clusters. The gas in galaxies is 
much denser, and the cooling times are, therefore, 
much shorter, and one might have expected that 
k might have a much more complicated profile. 
We demonstrate below that indeed the hot gas in 
galaxies does have virtually universal temperature 
ratios, with icfsl, 

3. Numerical simulations 

The hydrodynamical simulations were per- 
formed with GASOLINE (Wadsley et al. 2004; see 
Governato et al. (2007) or Stinson et al. (2010) 
for a more detailed description) — multi-stepping, 
parallel TreeSPH TV-body code. We include radia- 
tive and Compton cooling for a primordial mixture 
of hydrogen and helium. The star formation algo- 
rithm is based on a Jeans instability criteria (Katz 
1992), where gas particles in dense, unstable re- 
gions and in convergent flows spawn star particles 



at a rate proportional to the local dynamical time 
(see Governato et al. 2004). The star forma- 
tion efficiency was set to 0.05, but in the adopted 
scheme its precise value has only a minor effect 
on the star formation rate (Katz 1992). The code 
also includes supernova feedback (Stinton et al. 
2006), and a UV background (Haardt & Madau 
1996; Heller 1999). Additional simulations using 
the FTM 4.5 code (Heller & Shlosman 1994) are 
described in Section 3.1. 




bg(Density [cm 3 ]) 

Fig. 1. — Contour plot of the gas density and tem- 
perature of the individual gas particle at z = 
from the numerical simulation using GASOLINE of 
the galaxy G2. The gas particles either cool down 
to the floor of the atomic cooling 10 4 K, or stay 
at the DM temperature around 10 6 K. 

For the simulations with gasoline, we selected 
three candidate haloes with masses similar to the 
mass of the Milky Way (M « 10 12 M Q ) from an 
existing low-resolution DM simulation (300 3 par- 
ticles within 90 Mpc) and re-simulated them at 
higher resolution. These high-resolution runs are 
8 3 times more resolved in mass than the initial 
ones and included a gaseous component within 
the entire high-resolution region. In these runs 
the masses of the DM and gas particles are re- 
spectively nid = 1.17 x 10 6 /i _1 M Q and m g — 
2.3 x 10 5 /j _1 M Q . The DM has a spline gravita- 
tional softening length of 500 h^ 1 pc and there are 
about 10 6 particles for each component (dark and 
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gas) in the high-resolution region. More properties 
of the galaxies are listed in Table [TJ These hydro- 
dynamical simulations are described in details in 
Schewtschenko & Maccio (2010). 

Table 1: Galaxies parameters 
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Fig. 2. — Histogram of the individual gas particle 
temperatures at z = from the numerical simula- 
tion using gasoline. The three curves represent 
each of the three galaxies simulated at high reso- 
lution. The gas particles either cool down to the 
floor of the atomic cooling 10 4 K, or stay at the 
DM temperature around 10 6 K. 

In Fig. [1] we show contours over the temper- 
ature and local density of all the gas particles 
in one of the 3 galaxies (the others look essen- 
tially identical). In Fig. [2] we show histograms 
over the individual gas particle temperatures in 
the three selected galaxies, which is a projection 
of figure [T] The gas temperature exhibits a bi- 
modal distribution: it either cools down to the 
floor of atomic cooling about 10 4 K in dense clumps 
(central disk/bulge or satellite cores) or stays at 



the DM temperature in the larger and less dense 
regions. This is largely caused by the signifi- 
cantly short e r cool ing time around 10 5 K (see e.g. 



Keresetal 



(120051^ . When we discuss the gas 
temperature from the simulations using gasoline 
we refer to the second (large) bump. Practically 
we calculate a mass-averaged temperature by in- 
cluding only particles warmer than 10 4 5 K. Mak- 
ing a mass-average over all particles (including the 
cold ones) makes a negligible difference at most 
times. There is no metal cooling in this simula- 
tion, which is the reason for the floor at 10 4 K. 
The galaxies are typical field galaxies, and have 
been selected to have a quiet merger history (ex- 
cept Gl, which had a major merger at z = 0.8). 




Fig. 3. — The temperature ratios as a function 
of radius for the three simulated galaxies. The 
frames show that outside the central region of 
~ 6 kpc where the stellar disk dominates the 



potential, the ratio, n(r) 
proaches unity. 
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We can now consider radial bins for each galaxy 
at z = 0, and in each bin we calculate the mass- 
averaged gas temperature and the mass-averaged 
DM velocity dispersion. Their ratios as in eq.(|2]) 
is shown in Fig. [3] From this figure we see that 
outside the range of the disk, i.e. beyond ~ 6 kpc, 
these ratios remain close to unity. Beyond 20 kpc, 
where the gas density is lower and hence the cool- 
ing time longer, the gas temperature is slightly 



4 



above the DM dispersion velocities by 10%-20%. 
The velocity anisotropy is nearly constant, j3(r) ~ 
0.2 outside 2 kpc, slowly increasing to about 0.3 
near 100 kpc. We thus confirm that indeed there 
remains a substantial component of the hot gas at 
the DM temperatures even in dense environments 
of galaxies. 




0.4 0.6 0.8 

expansion factor a 

Fig. 4. — History of 2 hot gas particles. Shown is 
entropy, local density, and temperature as function 
of expansion factor. 

3.1. History of hot and cold particles 

In order to better understand the origin of the 
above bimodal distribution of gas temperatures, 
we tracked the evolution of various properties of 
individual particles. Each gas particle has its own 
unique thermal history, however, it appears that 
we as a first approximation can split the particles 
into four groups. 

• 1) Hot particles, which are shock heated to 
roughly 10 6 K, and then keep that temper- 
ature. 

• 2) Warm particles, which are shock heated 
to roughly 10 6 K, and then slowly cooled 
down somewhat. These constitute less than 
15% of the particles with T > 10 4 5 K. 

• 3) Cold particles close to the centre (r < 15 
kpc). These are shock heated, and then they 



cooled down very quickly. These may have 
experience several episodes of heating and 
cooling. 

• 4) Cold particles which are never shock 
heated. These particles come through the 
cold accretion phase. 

The two cold modes (3 and 4 ab o ve) were dis- 
cussed in detail in lKeres et al.l (|2005l) ; iMaccio et al 
(120061) . 

In Figs. (I4I5[) we present the entropy, local den- 
sity near the particle, and temperature as function 
of the expansion factor. In Fig. @] we present the 
history of 2 particles which are hot today. The 
first particle (Hot 1) only has one period of heat- 
ing as it enters the final halo. The second gas par- 
ticle (Hot 2) has one episode of heating followed 
by rapid cooling in its first subhalo, followed by 
heating when it enters the final halo. In Fig. [5] we 
present the history of two particle which are cold 
today. The first (Cold 1) has been accreted cold, 
and it never experienced any shock heating. The 
second particle (Cold 2) had 2 periods of heating 
followed by rapid cooling, which is clearly seen to 
happen when the particle enters the high density 
region. 




0.3 0.4 0.6 0.8 1 

expansion factor a 

Fig. 5. — History of 2 cold gas particles. Shown is 
entropy, local density, and temperature as function 
of expansion factor. 
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3.2. Evolution of k with redshift 

In order to test the generality of the above 
results, that n{r) ~ 1, we present results from 
i ndepen d ent si mulations of iRomano-Diaz et al. 
(|2008al lb1. 120091) using the parallel version of the 
hybrid SPH/jV-body FTM 4.5 code (Heller & 
Shlosman 1994; Heller, Shlosman & Athanas- 
soula 2007). The gravitational forces are com- 
puted using the falcON routine (Dehnen 2002) 
which scales as O(N). The tolerance parameter 9 
is fixed at 0.55. The gravitational softening ap- 
plied is e = 500 pc for the DM, stars, and gas. 
The vacuum boundary conditions are used and 
the simulations are performed with physical coor- 
dinates. The cosmological constant is introduced 
by an explicit term in the acceleration equation. 
The conservation of the total angular momentum 
and energy within the computational sphere in 
the colllisionless models is within ~ 0.01% and 
~ 1%, respectively. The evolution of various pa- 
rameters characterizing the DM and baryons is 
followed in 1000 snapshots, linearly spaced in the 
cosmological expansion parameter a. 

The modeling of star formation (SF) processes 
and associated feedback are described in Heller & 
Shlosman (1994) and Heller et al. (2007), which 
should be consulted for details. Feedback from OB 
stellar winds and supernovae (SN) Type II is also 
included. A fraction of this energy is thermalized 
and deposited in the gas in the form of thermal 
energy, then converted to kinetic energy through 
the equations of motion. This method is prefer- 
able over injecting a fraction of the stellar energy 
directly in the form of a kinetic energy (Heller et 
al. 2007). 

Multiple generations of stars are allowed to 
form from each gas particle. The evolution of gas 
metallicity is followed and the fraction of massive 
stars that lead to the OB stellar winds and SN is 
calculated from the Salpeter IMF. 

Fig. [6] displays the redshift evolution of the ra- 
dial profiles of the temperature ratio, k. All the 
gas particles (both hot and cold) are included in 
the averages here, and therefore different aspects 
of the time evolution of n(r) can be followed. At 
higher redshifts, z > 2, when the disk is being 
assembled from the cold gas supplied by the fila- 
ments, and which is converted into stars, the tem- 
perature ratio drops below unity in the disk re- 
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Fig. 6. — Evolution of n(r) with redshift using 
cosmological simulations of galaxy formation and 
evolution by Romano-Diaz et al. (2009). From 
top to bottom the frames show n(r) at redshifts 
z = 2,1,0. 

gion. For z < 1, when the star formation, feed- 
back by SN and stellar winds (Romano-Diaz et al. 
2008b, 2009) consume a large fraction of this gas, 
the temperature ratio rises to unity, as the heated 
gas tends to the DM temperature. After z « 1 we 
observe an influx of subhalos along the filaments 
which ablate the cold gas from the disk nearly 
completely. Two independent parameters verify 
the overall evolution of the model from a late-type 
to the early-type: the fraction of the gas-to-disk 
mass ratio, and the ratio of a spheroidal-to-disk 
mass ratio (within the disk radius). At z = we 
find a temperature ratio which is slightly below 
unity in the central region, and rises to slightly 
above unity towards the virial radius. 

Hence, we see that the gas-to-DM temperature 
ratio approaches unity for the evolved and relaxed 
galaxies. Of course, the actual evolution of the 
ratio does depend on the particular history of the 
galaxy, but the important point here is, that for 
equilibrated galaxies today, the temperature ratio 
is close to unity. 

As a sidenote we point out that the stellar dis- 
persion does not agree with the DM dispersion for 
these galaxies. We can therefore not use the stellar 



G 



velocities in the same way to extract information 
about the DM, and we cannot assume that the 
stellar velocity dispersion anisotropy should equal 
that of the DM. 

4. AGN outflows 

The simple relation between the gas and DM 
velocity dispersions in eq. (|3|) has been demon- 
strated us ing numeri c al sim ulations for galaxy 



clusters in Host et al 



(2009). Those simulations 
had sufficient resolution to probe a region outside 
of roughly 1 00-200 kpc. It remains unknown if this 
simple relation still holds further towards the cen- 
tral region, where cooling may be faster, or where 
a central AGN may pump energy into the intra- 
cluster gas. 

To address this question, we present results 
of a simulation of a three-dimensional model 
of AGN self-regulation in a cool-core cluster 
( Briiggen fc Scannapieco 2009). This simula- 
tion was per f ormed with FLASH version 3.0 



(|Frvxell et al.l |2000|) a multidimensional adap- 



tive mesh refinement hydrodynamics code, which 
solves the Riemann problem on a Cartesian grid 
using a directionally-split Piecewise-Parabolic 
Method (PPM) solver. In this simulation the 
cluster properties, such as the density and tem- 
perature profiles have been chosen to resemble 
the best-studied cluster, the Perseus cluster. The 
gravitational potential is taken to be static, and 
initially the gas is set up in hydrostatic equilib- 
rium. The gas physics includes radiative cool- 
ing, heating through AGN feedback and a subgrid 
model for Rayleigh- Taylor-driven tu r bulen ce, as 
described in Briiggen fe Scannapieco! ( 20091 ). Un- 
like in previous simulations, the energy input rate 
was not predetermined, but instead calculated 
from the instantaneous conditions near the centre 
of the cluster. A fixed fraction of the mc 2 rest 
mass energy of the accreted gas is returned to the 
ICM in the form of AGN-blown cavities whose 
turbulent evolution couples th em to the inflowing 
cool gas ( Briiggen et al. 20091) . When turbulence 
is properly accounted for, the model results in a 
sclf-rcgulatcd AGN, where a phase of an episodic 
heating is followed by a quiescent phase, where 
the cool-core reforms. 

For a wide range of feedback efficiencies, the 
cluster regulates itself for at least several 10 9 years. 



Heating balances cooling through a string of out- 
bursts with typical recurrence times of around 80 
Myrs, a timescale that depends on global cluster 
properties. 




radius [kpc] 



Fig. 7. — The temperature ratio as function of 
radius for the simulated galaxy cluster with a re- 
current AGN heating source. The full line shows 
the temperature ratio of the fiducial run of a 
self-regulating AGN in a cool core cluster after 
1 Gyr , labeled D5-10 in Briiggen fc Scannapieco! 
(2009). This figure demonstrates that outside 
the central region (« 30 kpc, where the recent 
AGN burst has heated the gas temperature much 
above the DM temperature) the ratio n(r) = 
<7 baryons/( CT DM) i s umt y, within 10%. The dashed 
line shows the results from an extreme case of 
AGN-fee dback, where strong outflows o ccur every 
50 Myrs (jScannapieco fc Bruggenll2008l ). 



In Fig. ([7]) we see that for the self-regulating 
AGN feedback model after 1 Gyr (solid line), 
the spherically averaged gas dispersion does in- 
deed agree with the DM dispersion within 10%, 
at radii larger than 30 kpc. This case rep- 
resents the best available simulation of a self- 
regulating cool core cluster which reproduces a 
number of key features seen in X-ray observations 
( Briiggen fc Scannapieco 2009t ). 

In the central region, at radii smaller than 
30 kpc, we see a sharp increase in the gas tem- 
perature resulting from a recent AGN burst, and 
at the very centre, within less than 10 kpc, we ob- 
serve the effect of very rapid gas cooling, and the 
re-emergence of a cool core. 

To test an extreme case of AGN heating of a 
cluster core, we also present the results from a 
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model where AGN-driven bubbles occur at reg- 
ular intervals of 50 Myrs with bubble energie s 
of 5 x 10 59 ergs (jScannapieco fe Briiggenl 120081 ) . 
Even for this non-self regulated, and somewhat 
less physically motivated case, the dashed line in 
Fig. ([7]) shows that, as long as we consider radii 
outside of 70-100 kpc region, the velocity disper- 
sion ratio stays within 10% of unity. 

5. Implications for the velocity anisotropy 

For the self-regulated AGN, the simple temper- 
ature relation k(t) ~ 1 holds outside the inner 
region of r ~ 30 kpc, as we have shown in sec- 
tion 4. For the o bserved X-ray clusters studied in 
Host et al. I <l2009h . this radius corresponds roughly 
to the centre of the innermost bin. This shows 
that one can trust the inference of the DM veloc- 
ity anisotropy in galaxy clusters down to radii as 
small as » 50 kpc. 

In the case of galaxies, we have demonstrated 
that the simple temperature relation holds for re- 
laxed and equilibrated structures, when applied to 
radii beyond the regions dominated by the stellar 
component. This opens for the possibility of indi- 
rectly determining the DM velocity anisotropy in 
galaxies in the following way. By measuring the 
temperature and the density of the hot X-ray emit- 
ting gas, one may use the equation of hydrostatic 
equilibrium, eq. (JTJ, to measure the total gravitat- 
ing mass. Now, considering the Jeans equation, 
eq. ©, we see that it contains the total mass on 
the left hand side, and 3 unknown quantities on 
the right hand side (pdm, of, ft)- First, combin- 
ing the total mass measurement from X-ray with 
spectral determination of the stellar mass, allows 
one to determine the DM density profile 

Pdm{t) = fi (T g (r),p g (r) 7 p stal {r)) . (5) 

We are thus left with 2 unknowns on the right 
hand side of the Jeans equation, eq. ([2]), namely 
the radial dispersion and the velocity anisotropy, 
(3, Using the temperature relation, eq. ([3]), to- 
gether with the X-ray determined gas tempera- 
ture, we get rid of one of these, and hence get an 
indirect measurement of the dark matter velocity 
anistropy 



0) = h (K(r),T g (r), p g (r), p star (r)) 



(6) 



ter (Hansen &: Piffarettal2007t iMorandi fc Ettori 
20071 iHost et al.l 12009^ and it is clear that such 
indirect observations in galaxies most likely will re- 
quire improved X-ray satelites. The simplest pos- 
sibility would be if re(r) = 1 everywhere, however, 
the method to infer the DM velocity anisotropy 
as described above, w orks as long as t he radial 
form of n{r) is known ( Host et al.l 12009). This is 
exactly one of the points of this paper, that the 
form of k(t) appears to be universal for galaxies 
(see figure [3]). 

A direct measurement of this velocity anisotropy 
is expected t o be virtually impossib le in terrestrial 
experiments ( Host fc Hansenll2007l ). Measuring /3 



of the DM in galaxies, as described above, can 
provide an independent confirmation to virtually 
every numerical modeling of DM halo formation 
in the cosmological context, and to the suggestion 
that a non-zero j3 is a fundamental pro perty of 
galac t ic DM halos, even in equilibrium ([Hansen 
120091 lHansen et all l2010h . Furthermore, a non- 
zero /3 has an effect o n the underground direct de- 
tection experiments (Vergados l2000l : lEvans et al 



12000 : IVergados et al 



2008). A detection through 



the method laid out above would therefore de- 
crease the systematic error-bars in terrestrial DM 
detection experiments. 

6. Conclusions 

In summary, we have found a very simple re- 
lation between the temperature of the hot gas 
and the averaged dispersion of the dark matter, 
namely that their ratio is close to unity, k « 1. 
We have demonstrated that this near equality of 
the gas and DM temperatures holds in the case of 
equilibrated and relaxed galaxies. 

This relation is not only conceptually impor- 
tant, but it will also allow an indirect determi- 
nation of the dark matter velocity anisotropy in 
galaxies. Such an observation will provide an im- 
portant and independent confirmation of all cos- 
mological DM halo formation simulations. 

We have also studied this relation in galaxy 
clusters containing an active AGN, and our results 
confirm that one can indeed use galaxy clusters to 
infer the dark matter velocity anisotropy. 



Such observation has only very recently be- 
come possible for X-ray bright galaxy clus- 
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